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Abstract 

We present a new next-to-leading order calculation for fully differential single-top- 
quark final states. The calculation is performed using phase space slicing and dipole 
subtraction methods. The results of the methods are found to be in agreement. 
The dipole subtraction method calculation retains the full spin dependence of the 
final state particles. We show a few numerical results to illustrate the utility and 
consistency of the resulting computer implementations. 



1 Introduction 



Single-top-quark production provides an excellent opportunity to study the charged- current 
weak-interaction of the top quark |]T]-§]. Measurement of the production cross section of 
single top quarks is planned at both the Fermilab Tevatron and the CERN Large Hadron 
Collider (LHC) |J. The DO |TIJ and CDF [II] collaborations have already set limits on 
both the s-channel and t-channel cross sections using data collected during run I of the 
Fermilab Tevatron, and have developed strategies for discovery at the current run II. 

Within the standard model, the single-top-quark cross section allows a direct measure- 
ment of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element Vtb- Further, the short 
lifetime of the top quark presents an opportunity to observe the polarization of the top 
quark at production, and hence directly probe the V — A nature of the weak interac- 
tion|TJ,|rj. 

Further, the mass of the top quark is comparable to the electroweak symmetry breaking 
scale, which has lead to speculation regarding new physics involving the top quark for some 
time ||14|| . There are many proposals to use these measurements to study non-standard 
couplings fi~5H23||, strong dynamics |24j-|27|l, flavor changing neutral currents |28|, CP 
violation p0|-|32][ , supersymmetry |33|-|35|], i?-parity- violating supersymmetry |36|-p^[, and 



Kaluza-Klein modes 

The characteristics of single-top-quark production form a new class of benchmarks for 
testing Quantum Chromodynamics (QCD). Experimental comparisons to the calculated 
kinematic distributions can provide a handle on inputs to the calculation. For instance, 
this process involves directly the bottom quark parton distribution function. Currently 
this function is constrained only indirectly though global fits [H-I3| to data. The bottom- 
quark density is then calculated from the light parton densities [24-50] in the context of 
QCD. Including a process that depends directly on the bottom density into the global 
analysis would be of great value to LHC observables, as their typical scale will be much 
larger than the bottom mass where the evolution begins. 

Finally, single-top-quark production is a background to all signals with VF+jet or W+b 
as backgrounds. This background is significant in a number of Higgs search channels (for a 
and other new physics, such as i?-parity conserving |^2[ and violating 
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supersymmetry searches. 

In this paper we present and discuss new calculations of single-top-quark production 
at next-to-leading order (NLO) in QCD. Earlier calculations of the NLO single-top-quark- 
production cross sections exist in the literature. The first calculation [55] was for the double 
differential cross section and used small masses for the gluons and quarks to regularize 
infrared and collinear divergences. Mass factorization was performed in the Deep Inelastic 
Scattering (DIS) scheme. Subsequent NLO calculations for the s-channel [^,0 and t- 



channel p8|,|58| modes used dimensional regularization and expressed the semi-inclusive 
cross section in terms of the Modified Minimal Subtraction (MS) factorization scheme. 
The value of our work is that the results are fully differential (meaning experimental cuts 
and jet finders can be applied), the results contain spin information, and the results use 
standard methods and schemes. The calculational methods used are an instructive step 
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toward the computation of NLO corrections to the tbj production channel that should also 
be considered when studying single-top-quark production. 

The NLO cross section receives contributions from virtual corrections and real emission 
diagrams. Taken separately, both parts are divergent and therefore cannot be evaluated in 
a straightforward way numerically on a computer. Only the sum of the virtual corrections 
and the real emission contributions is finite after mass factorization. Writing a general- 
purpose NLO Monte-Carlo based program therefore requires the analytic cancellation of 
singularities before any numerical integration can be performed. 

The two main general methods to handle the cancellation of singularities without loss 
of information are the phase space slicing |59|-p6| and the subtraction p7-73| methods. In 
this paper we implement the phase space slicing method of one [B5| and two |56| cutoffs, 
and the subtraction-based dipole formalism of Ref. |[73|| . We find that the results of all 
methods agree. The dipole calculation uses helicity amplitudes and therefore contains the 
complete spin correlations of the participating partons. 

We organize this paper as follows. In the next section we present an overview of 
the amplitudes entering the calculation. Section [3] discusses the framework of the phase 
space slicing method and presents the attendant analytic results. Section [| discusses the 
cancellation of divergences within the context of the massive dipole subtraction method, 
and gives the results for all relevant amplitudes. This Section also contains a detailed 
discussion of scheme independence for different ways of handling 75. The analytic results 
of Sees. |] and f| are presented using different, but self-consistent, notations appropriate to 
their methods of handling of infrared divergences. Brief numerical results are presented 
in Sec. |^, along with a comparison of the methods. Our conclusions are followed by an 
Appendix containing the relevant one-loop scalar integrals. 



2 Overview 

The lowest-order Feynman diagrams are shown in Fig. |T]. They are commonly distinguished 
by the sign of the W boson momentum squared. The t-channel flavor excitation process 

u + b-^t + d, (1) 

occurs via the exchange of a virtual space-like P^-boson, and the Drell-Yan-like s-channel 
process 

u + d^t + b, (2) 

occurs via a virtual time-like P^-boson. In reaction ([[]) it is understood that we may 
replace the (u, d)-quark pair by (d,u), (c,s) and (s, c). In reaction (Q) we may replace 
the (u, d)-pair by (c, s). In addition, CKM suppressed combinations are included at all 
vertices. 

For each channel the fully differential spin-averaged Born cross section can be written 

as 

^ 0) = ^El^Bo r n| 2 rfPS 2 , (3) 
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(a) 



(b) 



Figure 1: The leading order Feynman diagrams for reactions given in (a) Eq. | and (b) Eq. §. 
The top quark line is doubled. 



where s is the partonic center-of-momentum energy squared, and we use p = s, t to denote 
the channel. The two body phase space is given by 

1 d 3 pi d 3 p 2 



dPSo 



'-S (4) (q -pi -p 2 ) • 



(4) 



(2tt) 2 2Ei 2E 2 

The t-channel Born matrix element squared summed (averaged) over final (initial) state 
spin and color states is 

^ |-Mt,Born| 



(5) 



Here, s = (p u + pb) 2 and t = (p u — Pd) 2 , the partonic reaction sub-energy squared and the 
square of the momentum transfer across the W, respectively. The CKM matrix elements 
\Vij\ may be changed for the given particles, and mt is the top-quark mass. The result 
for the s-channel is obtained by interchanging s and t, and letting s = (p u + pg) 2 and 
t = {Pu — Pi) 2 - I n Sec. |] we discuss the cross section without the sum over spins. 




(b) 

Figure 2: The one- loop virtual corrections to the (a) t-channel and (b) s-channel amplitudes. 



At next-to-leading order we must include the virtual QCD corrections to Eqs. (HD and 
, shown in Fig. 0. To obtain the real emission contributions one must attach a gluon in 



4 



all possible ways to the diagrams in Fig. [L| The resulting crossed diagrams may be written 

as 



u + d -> t + b + g , (6) 
u + b — > t + d + g 

d + b -> t + u + g , (7) 

d + g -y t + b + u, (8) 
5f + fe^t + rf + M. (9) 

These crossings may belong either to the s or t channel, and their assignment is governed 
by the collinear singularity structure. It is useful to distinguish whether the gluon attaches 
to the fermion current containing the heavy-quark line (h) in Fig. |] or to the current 
containing only the light-quark lines (I). Note that the contributions from h and I radiation 
do not interfere in the cross section due to color conservation 58|. 



The diagrams in Eq. (|5]) belong fully to the s-channel, and diagrams in Eq. (|7]) belong 
to the t-channel. For crossing in Eq. (|8|) let us consider the first process u + g — > t + b + d. 
In the heavy-quark line the gluon can split into a collinear 66-pair with reduced process 
u + b — > t + d, so that this crossing belongs to the t-channel. In the light-quark line the 
gluon can split into a collinear dd-pair, with reduced process u + d — > t + b, so that this 
crossing belongs to the s-channel. The classification of the second process in Eq. (|8|) is 
similar, with the role of the u- and <i-type quarks exchanged. Finally, in Eq. @, the gluon 
connected to the light-quark line can split either into a collinear wM-pair or into a collinear 
dd-pair. In both cases the process reduces to the t-channel process u + b — > t + d or 
d + b — > t + u. Crossing of Eq. (^|) to the heavy-quark line is not included in our calculation 
as it involves on-shell H^-decay into light fermions, and is classified as VF-top associated 
production. The cross section for this process is estimated to be negligibly small for the 
Tevatron 0,§,|75|-|77|]. 



In what follows we consider only the production of a single top quark, but the relevant 
expressions for the production of an antitop quark may easily be obtained by charge con- 
jugation. Our calculational framework is further specified as follows: we set the bottom 
quark mass to zero, so that we work with 5 massless flavors in the parton distribution func- 
tions. It is important to use a consistent set of 5-flavor parton distribution functions, in 
which the 5-flavor set has been computed from a lower flavor number set via NLO matching 



conditions |^|| that preserve the momentum sum rule. To handle divergences occurring 
at intermediate stages of the calculation we use dimensional regularization (D = 4 — 2e) . 
Details describing renormalization of the vertices that contain 75 are given in Sec. [4.2| . 
Collinear divergences are subtracted in the standard MS scheme. We use the Feynman 
gauge for the gluon propagator. 
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3 Phase space slicing 



In the phase space slicing method a sub-region of phase space containing soft and collinear 
singularities is defined. If the subregion is small enough, one may make simplifying kine- 
matic approximations in the phase space integrals, so that they can be performed analyti- 
cally. To define the size and shape of this subregion one introduces one or two theoretical 
cutoff parameters (for common or separate treatment of collinear and soft contributions, 
respectively). The final cross section should not depend on the choice made for the parti- 
tioning. Below we discuss the application of these methods to single-top-quark production. 



3.1 Phase space slicing with two cutoffs (5 S ,5 C ) 

A detailed description of the phase space slicing method with two cutoffs has recently 
been given in Ref. ||66| . In this section we present the results needed for single-top-quark 
production. We follow closely the notation of Ref. . 

In the two cutoff method, phase space is divided into a hard (H) and soft (S) region. 
The contribution from the latter region is computed by applying the eikonal approximation 
to the radiative processes in Eqs. (§]-fJ). The soft region of phase space is defined by a 
condition on the energy of the soft gluon in the partonic center-of-momentum frame: 



< E g < 6 S ^ 
2 



(10) 



The hard region is then defined by E g > 5 s ^/2. After performing the D-dimensional 
angular integrations of the gluon in the region defined in Eq. (|10"D , the soft contribution is 
found to be 



daf = da® 



a s T(l - e) /4tt/z 2 



27rr(l - 2e) V s J 
with p = s,t labeling the channel. For the s-channel 



A\ 

An 



+ 



3C f 



c h 



6 In 5, 



In ( — 2 



6 In 2 5 S - 2 ln<L + 2 ln<Lln — 



s + m 



In 



77T 



rrr 



-2Li 2 (/3)--ln 2 



rrr 



'121 



and for the t-channel 
3C F 



A\ 



A\ 



6 In 5. 



In 



m 2 s 



6 



An 



c h 



6 In 2 8 S - 2\n8 s + 41n5 s ln 



s(3 



+2 In S s In 



[m 



m 2 s 



+ 



s + m 



In 



s — 



+ ln 2 



s(3 



2Li 2 1 



s3 



m 



m 2 - t 



2Li 2 — 



2Li 2 



s + u 



(13) 



where the top-quark mass is denoted as m, and (3 = 1 — m 2 / s. 

The hard region is further divided into hard collinear (HC) and hard-noncollinear (HC) 
regions. The latter is computed in 4 dimensions, integrating numerically over the HC phase 
space using standard Monte Carlo methods. The HC contributions arise from integration 
over those regions of phase space where any invariant, = (pj + Pj) 2 or = (pi — Pj) 2 , 
appearing in the denominator becomes smaller in magnitude than 5 c s. The singular regions 
are distinguished by whether they come from initial or final state radiation. The latter are 
given by 



with 



A 



47T/T 



C F [2\n6 s + 



2nT(l - 2e) 
3 



e 



(14) 



Z-^-ln 2 , 
2 3 



2 In/? 

\n 2 (3 + 2\n5 s ln/3 



\n6 c (2ln5 s + 



2 ln/3 



(15) 



for both channels (all massless partons in the Born cross sections are fermions). 

The collinear contributions from the initial state are given by the sum of two contribu- 
tions. The first is the finite remainder after mass factorization. The second results from a 
mismatch in integration limits when subtracting the mass factorization counter-term. 



7 ij->FS 
d(T p,C 



da® 



T(l-e) 



'47T/i 2 



2vrr(l - 2e) 



f?{ z > Vf) + ( — + A S q ) ff (z, fJ, F 



cH 



(16) 



The modified parton distribution function / is given in Ref. [66]. For the reaction at hand 
we only need consider the quark-quark initial state splitting, so 

A{ c = C F (2ln5 s + 3 



ASC 

^0 



C F [2\n5 s + -) In ^ 



(17) 



The virtual contribution is obtained as explained in the Sees. |4.1| and f|72|. The results in 
the notation of this section are 



daW = daf 



) /47T/X 2 

27rr(l -2e) I s 



'AX 



AV n 



^ ) da™ 
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where 



A\ 



Cf{[-2] 

C p 

c F 



[1]} 

2 In 



+ 



5 
2 



2 In (1 - A) - In 



In 2 



3 In 



7T 

IF 



+ 



m 



^ln( — 

m 



2 In (1- A) In — 



A 



In (1 - A) - In 2 (1 - A) - 2 In (1 - A) + 2Li 2 (A) 



-6 

TT 2 

y 



(19) 



In the above we have defined A = q 2 / (q 2 — m 2 ). Further, the separate terms in the square 
brackets originate from the massless-massless or the massive-massless vertex corrections. 
For the t-channel one sets q 2 = t, while for the s-channel one sets q 2 = s. 

, 2 



1 1 



g \V U d\ \V t b\ C F 



2n m 2 _su^ 
t 



m 



1 



m 2 — t 



(20) 



is the t-channel finite piece in the virtual contribution that is not proportional to the Born 
cross section. It results from the interference of the renormalized massive-massless vertex 
with the Born amplitude. The s-channel version may be obtained by crossing. 

At this point one can see that the two body weight is finite: A 2 + A\ = and A\ + 
A\ + A\ + 2A\ C = 0. The factor of two occurs since there are two quark legs, either of 
which can emit a gluon. The final finite two-body cross section is given by the sum of the 
residual / terms from both the quark-quark and quark-gluon initiated processes and the 
finite two-body weights. The result, summed over parton flavors is 



a 



(2) 



2^ 



J dx 1 dx 2 {fa 1 (x 1 ,n F )f b -{x 2l HF 

a,b 

+ H 0) [f^ 1 (x 1 ,ix F )f b H2 (x 2 ,ix F ) + f? 1 (x uf x F )f* I *(x 2 ,v F ) 
The three-body contribution is given by 



daf ] (A p + A% + A + 2A 
x 2 )) 



(21) 



a 



(3) 



with 



a, 6 



(aft) | 2 



J2ldx 1 dx 2 -J H _Y,\M 3 



(aft)|2 d p 



-7ig 4 \V ud \ 2 \V t 



tb\ 



3 • 



1-3 



(22) 



(23) 



The contain the luminosity and Dirac algebra of Eqs. ([T^D . We choose initial-state 
momenta as incoming, and label the momenta for the t-channel as 



#1 : u(pi)b(p 2 ) 
^2 : u(pi)g(p 2 ) 
^3 : g(pi)b(p 2 ) 



d{pz)t{p±)g(ps) , 

d(p3)t(p4)b(p5) , 
d{pz)t{p A )u{p 5 ) . 



(24) 
(25) 
(26) 
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The t-channel are given by 

s 12 (^14 — s u) — s 34^25 S 12 ( s 34 + s 45) — s 34^23 



*a = 2C F F l 



*15 s 35 



_ ^13 [ S 12 (2^34 + S45) + s' 34 t 25 

+2C F F 2 



^15 s 35 / 
s 12 (^23 — S 34) — s 34^15 S 12 ( S U + S 35) — 2m 2 /s 45 ) ~ S uA 
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^25 S 45 



^24 [ g 12 (2-^34 + S 35) + fgAj)] 
^25 S 45 



(27) 



, T , 77 I s 12 s 34 + ^15 ( S 34 s 35) , ^15 ( S 34 + ^23) (1 2m t /t' 24 ) S 34 t' 14 

^2 = ^2 " + 

\ r 25 r 24 

1 S 45 [ S 12 S 34 + ^15 (^23 + ^ S 3i)\ \ fr,o\ 

+ vi ' { ' 

^ - P ( g 12 S 34 + ^25 ( g 34 ~ ^45) , ^25 ( g 34 + ^14) ~ g 34^23 
1 S 35 [ s 12 s 34 + ^25 (^'l4 + 2s' 34 )] \ 

+ p: , (29) 

t 13 t 15 / 

where C F = 4/3, s^- = (p, + p.,) 2 , ^ = (p* - p^) 2 , s^- = s^- - m 2 , ^ = ^ - m 2 , 

Fi = i^kf u ■ (30) 
*-fe^sw*- (31) 

a s i(h) and the luminosity functions L l{h) = f^ Il (x 1 , fi Fl{h) )f^ 2 (x 2 , fi F i{h)) are evaluated 
using the scales at the light (heavy)-quark lines, respectively. All other s- and t-channel 
matrix elements can be obtained by crossing. Physical predictions follow from the sum 
+ a^ 3 \ which is cutoff independent for sufficiently small cutoffs as shown below. 



3.2 Phase space slicing with one cutoff (s m i n ) 

The calculation using the one cutoff slicing method is similar to the one using the two 
cutoff slicing method, with some differences that we now address. In this method, a pair 
of partons with momenta p« and Pj is defined to be unresolved if 

|2pi -Pj\ < S min , (32) 

with s min small compared to the hard scale of the process. The condition in Eq. (|3"2"D 
can occur if either pj and pj are collinear, or if one of the two is soft. This method, 
combined with the use of color-decomposed amplitudes and universal crossing functions, 
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has been developed into a general method for computing with minimal calculational effort 



fully differential NLO production cross sections of bosons and jets in [|63i|64| and identified 
hadrons and heavy quarks in JHJ . The single-top-quark production process has a relatively 
simple color structure, so we do not need to decompose the scattering amplitudes into 
color-ordered subamplitudes. 

The treatment of the virtual corrections is no different from the two cutoff method. 
To determine the radiative corrections, all partons are first crossed to the final state, and 
resolved and unresolved contributions are identified according to the criterion in Eq. ([32]). 
The unresolved contributions, (soft and collinear) can be found in [B3-65J expressed in 
D = 4 — 2e dimensions. The soft contributions are expressed in terms of the momenta for 
partons in lowest order kinematics with all partons in the final state 



-> u + d + t + b. 



and are given by 



a s C F 



'47T/i 2 



7T r(l-e) 



e 1 V s ■ 



+ J(m, 0) 



Pb 



(33) 



(34) 



with da^ obtained by crossing all momenta to the final state, and 



1) 2p t -p b >m 2 : J(m,0) 

2) 2p t -p b <m 2 : J(m,0) 



1 1 (2p t -p b 



2e 2 V 



J_ (2p t -p h 



2e V 



7T 



m 



( 2pfPb 
V m? 



+ 



7T" 



2e 2 2e 12 



+ 1 • 



12 2p t ■ p b 
(35) 



The s-channel contribution is then obtained by replacing 2p u ■ p^ —>■ s, 2p t ■ p b —* s — m 2 . 
The t-channel contribution is obtained by replacing 2p u ■ — > t, 2p t ■ p b — > t — m 2 , which 
leads to n 2 terms after expanding in e. The collinear contributions are likewise given by 



da^ 



-da® 



x 




(36) 



m 



|. The s-channel contribution is obtained by 
The t-channel contribution is obtained by 



The / functions are given in Refs. |5B 
replacing 2p u ■ p d -> s, 2p t ■ p b -> s - 
replacing 2p u ■ pa —* t, 2p t ■ p b —>■ t — m 2 . The sum of these contributions is already finite. 
One now generates the various subprocesses of single-top-quark production by crossing 
pairs of partons back to the initial state. Crossing symmetry is not a property of next-to- 
leading order cross sections, but it may be implemented in the following way. 

In general a NLO fully differential cross section for a process with initial hadrons H\ 
and H 2 may be written as 



da HlH2 = Yjf dx lf dx 2Fa l { X l)Fb 2 { X 2)d°ab LO { X li X 2 
a,b 



(37) 
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where a, b denote parton flavors and X\, x 2 are parton momentum fractions. The function 
da^ b LO is computed with all-partons-in-the-final-state matrix elements, in which partons a 
and b have simply been crossed to the initial state, i.e. in which their momenta p a and p b 
have been replaced by — p a and — p b (this function does include the 7r 2 terms resulting from 
this replacement in the one-loop virtual graphs). The functions T^{x) are modifications 
of the parton distribution functions f^(x,fj, F ) 

?a{x) = f?(x,fi F ) + a s C?(x^ F ) + 0(a 2 s ), (38) 



where C^(x,n F ) are finite, universal "crossing functions" |64| |. They implement the cross- 



ing property for the unresolved contributions, and are given by 



Af(*,^)ln l^j+B^(x^ F ) 



(39) 



The functions A^, functions for the proton are given in |54|. In the unresolved contri- 



bution one may simply cross pairs of partons without further modifications. 
The full NLO differential cross section can now be written as: 

dcr HlH2 = / dx 1 dx 2 fa 1 {x 1 ,fi F ) f b H2 {x 2 ,fi F ) da^ b LO (x 1 ,x 2 ) + a s (/i F ) 

(Cf 1 (x u LL F )f b H2 (x 2 , MF ) + (x u fx F )C b H > (x 2 , ^ F ))da^ b (x 1 , x 2 )] . (40) 



X l 



The unresolved contribution, now including the crossing functions, depends analytically 
on s m i n , but this s m i n dependence cancels against that of the resolved contribution. The 
results produced with this method agree with those of the previous section. In this paper 
we limit ourselves to some illustrative numerical studies, so that we only employ the dipole 
and two-cutoff slicing methods for numerical results. 

4 Massive Dipole Subtraction Calculation 

Within the dipole formalism the NLO cross section is rewritten as 
a NLO = J da R + J da v 

n+l n 

= J (da R -da A ) + J (da v + J daA . (41) 

n+l n \ 1 ) 

In the second line an approximation term da A has been added and subtracted. This 
is valid if all singularities occur in the final state. For initial state partons there are 
slight modifications. The approximation da A has to fulfill the requirement that da A is a 
proper approximation of do R such as to have the same point-wise singular behavior (in 
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D dimensions) as da R itself. Thus, do A acts as a local counter-term for da R and one can 
safely perform the limit e — > 0. This defines a cross-section contribution 



„NLO 
a {n+l} 



( da R -da A ) . (42) 

V e=0 e=0J v ' 

n+1 

da A is analytically integrable (in D dimensions) over the one-parton subspace leading to 
soft and collinear divergences. This gives the contribution 



NLO 
a {n] 



J ida v + J da A \ . (43) 
n V 1 / ^=n 



e=0 

The final structure of an NLO calculation is then 

„NLO JVLO | NLO f A A\ 

o = (r {n+1} + cr {n} . (44) 

Since both contributions on the right hand side of Eq. ([0]) are now finite, they can be 
evaluated with numerical methods. The (n + 1) matrix element is approximated by a sum 
of dipole terms 

do A ~ £ J2v ijik 

pairs i,j k^i,j 

= E E - 9n 1 - ( 1 » ••■> •••> fc. Tfc T 2 T ' J ^l 1 ' fa')' ^ •••) » 

pairs fc^i ,j Z "* ' Pj 1 ij 

(45) 

where the emitter parton is denoted by ij and the spectator by k. Here Tj denotes the color 
charge operator [|72| for parton i and is a matrix in the helicity space of the emitter 



with the correct soft and collinear behavior. |1, (ij), ■ k, ...) is a vector in color- and 
helicity space. By subtracting from the real emission part the fake contribution we obtain 

da R -da A = # n+1 ^|M( R ,...,p n+1 )r^ 1 (pi,...,p„ + i) 

E E V ij,k(P^ -,Pn+l)dn Ut (PU -,Pij, -,Pk, -,Pn+l) J • 
pairs i,j k^i,j / 



(46) 



Both da R and da A are integrated over the same (n + 1) parton phase space, but it should 
be noted that da R is proportional to # n + 1; whereas da A is proportional to 9^. Here 6^ ni 
denotes the jet-defining function for n-partons. 

The subtraction term can be integrated over the one-parton phase space to yield the 
term 

1® da B = j ' da A = J2 E / d(f>di P oleV ij>k . (47) 

i pairs i,j k^i,j 

12 



The universal factor I still contains color correlations, but does not depend on the unre- 
solved parton j. The term I g) da B lives on the phase space of the n-parton configuration 
and has the appropriate singularity structure to cancel the infrared divergences coming 
from the one-loop amplitude. Therefore, 

da v + I <g> da B (48) 

is infrared finite and can easily be integrated by Monte Carlo methods. The explicit 
forms of the dipole terms P^fc, together with the integrated counterparts, can be found in 
Ref. JF2] (the original massless case) and in Ref. [|73| (extension to massive fermions). 



4.1 Calculation of the amplitudes 

We have performed three different calculations of the loop amplitudes. One calculation 
was done using the standard approach in the 't Hooft-Veltman scheme. The second calcu- 
lation involved treating the 75 as anti-commuting in D dimensions, thereby retaining Ward 
identities for the charged current vertex. In the third one we calculated helicity amplitudes 
using a four-dimensional scheme. The effects of different prescriptions for 75 are discussed 



more extensively in Sec. 4.2. The results of the three calculations agree with each other, in 



the sense that they can be related to each other through process-independent finite renor- 



malizations. In addition, we find agreement with the earlier calculations of Refs. |78|.[79 
We present here the helicity amplitudes obtained with the third approach. They are more 
compact and contain the complete spin information. In the standard approach one just 
calculates the interference between the loop amplitude and the Born term and sums over 
all spins. 



We first list our conventions (for reviews of spinor helicity methods see e.g. |80|,|8T 



With spinor helicity methods we can express scattering amplitudes in terms of massless 

1 

2> 

u(p, ±) = v(p, T) = |p±) , u(p, ±) = v(p, T) = (p ± I ■ (49) 



Weyl spinors of helicity ± 1 



External fermion states are directly expressed in terms of these. Our convention is to take 
all particles outgoing. For example, an outgoing massless fermion with positive helicity is 
denoted by (p + | , while an outgoing massless anti- fermion with positive helicity is denoted 
by \p— ). The gluon polarization vectors, of helicity ±1, may be written as 

We use the customary short-hand notation: 

(ij) = (Pi - \Pj+), M = (Pi + \Pj~)- (51) 

In Eq. ([5(]) k is the gluon momentum and q an arbitrary light-like "reference momentum" . 
The dependence on the choice of q drops out in gauge-invariant amplitudes. We shall also 
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employ the abbreviations 



(i-\k + l\j-) = (ik}[kj) + {il}[lj], 

Sij...k = (Pi + Pj + ■■■+ Pkf , (52) 

with all momenta null-vectors. 

To treat the massive top quark within the framework of spinor helicity methods, we 



use the extension to massive fermions |82|-|85|. Even though helicity is not a conserved 
quantum number for a massive particle, a massive positive-energy spinor satisfying the 
Dirac equation has a two-fold degeneracy (labeled by a spin-component quantized along 
some axis). With slight abuse of notation we label these two states by "+" and "— ". Let 
p be a four- vector with p 2 = m 2 and p$ > 0, and let q be an arbitrary null vector with 
go > 0. We define 

1 1 

u(p, +) = —~= {p + m) |g-), v(p, +) = —= (p - m) g-}, 
y/2pq y/2pq 

1 1 

u{p,-) = -=={p + m) \q+), v(p, -) = -=={p-m) \q+). (53) 

For the conjugate spinors we have 

u(p, +) = -^=(Q -Kp' + 'm), v(p, +) = -7=(g - | - m) , 
V2pq y/2pq 

= —7=(q+ \ty + m), v(p, -) = -==(q+\($-m). (54) 
y/2pq y/2pq 

It is easy to check that for these spinors the Dirac equations, orthogonality, and complete- 
ness relations hold. The dependency on the arbitrary reference momentum q drops out in 
the final answer. 

Given two four-vectors p and q, the spinor product (pq) is calculated as follows: If 
p t > and q t > 0, 

(pq) = — (p±q+ - p+q±) ■ (55) 

Here the light-cone coordinates p + = pt+p z and p± = p x +ip y are used. For negative-energy 
four-vectors we have 

(PQ) = -<(-?)?), forp t <0, 

(pq) = ~(p(-q)), forg t <0. (56) 
The spinor product [pq] is related to (pq) by 

[pg] = sign(p t q t )(qp)* . (57) 
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We employ amplitudes with all partons outgoing, generating in an economical way the 
relevant scattering amplitudes by crossing. The lowest-order amplitude is shown in Fig. [l 

A Wb : -> t(p 8 ) + b(pi) + d{p 6 ) + u{p 7 ) , (58) 

where all momenta are outgoing. We use here the notation of Ref. [5B|, which explains the 
unusual labeling of the momenta with p 4 , p 6 , p 7 and p$. 

Each amplitude we decompose into gauge-invariant partial amplitudes. The color de- 
composition of the Born amplitude is 

Awb,born = ^84^67 A\y bborn . (59) 

It is convenient to factor out some common prefactors from the partial amplitude A^ b born 
and to write it as follows: 

a[1] e2V ud V tb 2i B Wb,born / Rn \ 

^■Wbborn — r> • 2 n 9 / V 0U / 

Wb ' born 2sm 2 6 w s 67 -m 2 w ^ _ |4 + 6 + 7|2-) 

Here we denote the reference momentum for the massive spinor by q = p 2 . If one is 
only interested in the spin- summed squared amplitude, one may choose any arbitrary null 
vector for p 2 , the choices p 2 = (1, 0, 0, 1) or p 2 = p& are examples. However, by keeping p 2 
unspecified, we keep the complete spin information, and our formulas become only slightly 
more lengthy. The non-vanishing Born amplitudes are 

B§ b , born (pt,P^,ptPs) = [47](6-|4 + 7|2->, 

B$ bM . n (pj,P6,Pt,Pt) = m(26}[74]. (61) 

We now turn our attention to the loop amplitudes. The color decomposition of the 
one-loop amplitude is given by 

A X T^ffl TO, A W i TO, rpa r A I 2 ] _i_ r T a T a A P] 

SiWb,loop — <J67 1 Si 1 ii^WbJoop ' 1 Gi 1 i7 u S4^-Wb,loop ' 1 67 1 8A /± Wb,loop 

iV 2 -l . .m N 2 -l 



x s aK I ~ L x s 4 [2] 

2 N °w 6 M A wb,ioop H 2N ° 67 ° M/i Wb,ioop 

1 /- - 1 ~ ~ N ,[3] 



+ 2 ^64^87 - J A wbloop . (62) 

Here we used the short-hand notation 

= T? sJi , (63) 

where a is the color index of the gluon, i 8 is the color index of the quark t(p$) and j'4 is 

the color index of the quark b(p4). A^ b ^ loop corresponds to loop corrections on the t-b line, 

[21 [31 
Awb,ioo P t° corrections on the u-d line, and A\y bloop to a gluon exchange between the two 

[31 

lines. Note that we do not have to calculate A\^ bloop . 

2 Re (5 u 6 67 A$ b y T2 7 Tg 4 A® b>loop = 0, (64) 
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because the color matrices are traceless. We write 



2 sin 2 9 W s 67 - m* w ^-(2 - |4 + 6 + 7|2-> 



• ,r "-"-' - 2 sin 2 ^ S67 - w? w ^-(2 - |4 + 6 + 7|2-> W ' 



For the helicity configuration pt,p G ,p 7 ,Pg we obtain 
2[47]{(6 - |4 + 7|2-)(4 - |6 + 7|4-)^ a) ( S67) m 2 ) 



+ 



(6- |4 + 7|2-> + -- 5467 (64) [42] 



2 (4- |6 + 7|4-)/ N 1 1 ' 2(4-|6 + 7|4 
-(s 67 (6- |4 + 7|2-) + S467 (64)[42])4 f,) (m 2 ) 



B { a) (sQ 7 ,m 2 



(4- |6 + 7|4- 

I 1 (64)[42] A (m^ s (64) [42] ( _ 2g) 1 

+ 2(4-|6 + 7|4-) Ao(m } " S467 (4 -|6 + 7|4-p (Sfi7 ' m } / 

-^Ufc,/oop(P4 iP% 1V7 )P8 ) = 

-2[47](6 - [4 + 7|2-> (* 67 cf ( Se7 ) + ^ c) (s 67 )) ■ (66) 
For the helicity configuration p% ,pg ,p 7 ,p§ we obtain 

-^U b.ioop (^4 5 % > > Pt ) = 



-2m [74] 



(62)(4-|6 + 7|4-)^ a) ( S67 ,m 2 ) 



+ 



Ys 67 (62) - i(67)[74](42)VS a) (s 6 7,m 2 ) 



(4-|6 + 7|4-)V 

' ( S67 (62) + (64)(4+|6 + 7|2+))4V 2 ) 



(4- |6 + 7|4-) 



1 (64)(4+|6 + 7|2+) (2-|6 + 7|4-)(46) ( _ 2e) 

+ 2 S467 (4 -|6 + 7|4-) Ao(m } + (4-|6 + 7|4-) C ° (S67 ' m } 



-2m[74](26) (s 67 cf (s 67 ) + \s 67 ) ) . (67) 



The expressions for the standard scalar integrals are collected in the Appendix. The 
ultraviolet (UV) renormalization is discussed in the next section. 
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Finally, we need the real emission amplitudes with one additional gluon. These are 
listed in Eqs. (|6|-||), and correspond to the process 

A Wg : -> t(p 8 ) + 6(p 4 ) + g(ps) + d(p 6 ) + u(pr) . (68) 
The color decomposition reads 

A W g,real = T 8A^7 A^- g real + ^^^Ws.reoi ■ ( 69 ) 



These amplitudes have been calculated in Ref. ||86|| . For convenience we list them here. We 
have 

[i] 9e 2 V: d Vt b {-i)2y/2 B$ gireal 



Awg,real 



2 sin 9 W s 67 - m 2 w ^_^ 2 _ |4 + 5 + 6 + 7|2-) 

,[2] ge 2 V: d V tb {-i)2V2 B§ greal 

w ^ " 2sin 2 ^ 5567 -m^ / (2 - | 4 + 5 + 6 + 7|2-)' 1 J 



As reference momentum for the massive spinor we have chosen q = p 2 . The non- vanishing 
amplitudes are 

Bwg,real (Pt 1P5 iP& iPl iP%) 

(6- |4 + 5 + 7|2-) / (6- |4 + 5|7-) [74] (6 - |4 + 7[5— ) \ 
(65) ^ (45) + S467 - ™ 2 /' 

Bw g ,real fat i Pt i P§ i P7 i P&) 

_m(26) /(6-|4 + 5|7-) [74] (6 - |4 + 7|5-)\ 
(45) + S467 - m 2 J' 

Bwg^real (Pa 5 ^5 J ^6 ' P"l ' ^8 ) 

= [54] (,l 4 '- m .) ((^ - 14 + 6 + 7|2->[47](76> 4- ^[24]<56) 

E>Wg,real (Pa 5 2>5 ) P% 1 Pi 1 P8 ) 

m [47] 



s 467 - m 2 [45] 



((25) (67) [74] + (56) (2 - |5 + 6 + 7|4— >) , (71) 



Bwg,real (Pa i Pt i P& i Pi i P% 



Bw g ,real {Pa i Pt i P& i Pi i P% 



(6 - |4 + 5 + 7|2-)(6 - |5 + 7|4-) 

WW) ' 
m(62)(6- |5 + 7|4-) 

<56><75> ' 



B m ( V + V - V - P + P -) - [74](2+ 1(4 + 5 + 6 + 7)(5 + 6)|7-) 

D Wg,real \Pa 5 P5 5 PQ ' Pi > ^8 I ~ jgyj jggj J 

S M fp + „- v - „+ - ^[47]<2- |5 + 6|7-> 

D Wg,real\PAiPhiP%ii ) l iV?,) ~ [57] [56] ' 
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The matrix element squared is given by 



\A W g,Real\ = ^N C (N 2 C -1) 



A m 

^Wg,Real 



+ 



A P] 

A Wg,Real 



(73) 



There are no interference terms between A^ g Real and A^ g Real . 



4.2 Conversion between schemes and scheme independence 

The one-loop amplitudes presented in the previous section have been calculated in a 
four-dimensional scheme. They differ from the corresponding amplitudes in the 't Hooft- 
Veltman scheme by finite terms. These finite terms are either of ultraviolet (UV) or infrared 
(IR) origin and result from expressions of the form e/e. To obtain the unique and correct 
result, one proceeds through the following steps: 

• The bare one-loop amplitudes may contain UV-divergences. With the help of a spe- 
cific regularization scheme these divergences are isolated and removed by renormaliza- 
tion. The explicit form of the renormalization depends on the chosen renormalization 
scheme. 

• The specific combination of regularization and renormalization scheme may break 
certain Ward identities. These Ward identities have to be restored through finite 
renormalizations. The required finite renormalizations are universal, i.e. they do 
not depend on the process under consideration. After this step all finite parts of 
UV-origin are uniquely fixed. 

• In addition, QCD amplitudes may contain IR-divergences. Unitarity requires that 
we employ the same regularization scheme in the phase space integral over the unre- 
solved real emission part as in the one-loop integral. 

Alternatively, since the structure of the IR-divergences is universal, we may derive 
simple formulas, which relate the finite parts specific to a certain regularization 
scheme to the ones of another scheme. We can therefore convert a one-loop ampli- 
tude calculated in one scheme to the corresponding amplitude calculated in another 
scheme. 

We will discuss the three steps in detail for the case at hand. To start, let us briefly 
summarize the properties of the 't Hooft-Veltman scheme and the four-dimensional scheme. 
The 't Hooft-Veltman scheme treats unobserved particles (particles in loops and unresolved 
partons in the real emission part) in D = 4 — 2e dimensions. Observed particles are taken 
in four dimensions. 75 is a four-dimensional object in the t' Hooft-Veltman scheme, anti- 
commuting with the first four Dirac matrices and commuting with the remaining ones. 

The four-dimensional scheme is specified in simple terms by the fermion propagators 
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Four- dimensional Dirac-matrices occur in the numerator, whereas D-dimensional quantities 
occur in the denominator. Two adjacent Dirac-matrices in the numerator are contracted 

as 

i>{4)i>{A) = (v\ D ) -V(-2e)) ■ 1(4) , (75) 

which can be interpreted as the statement U D is effectively larger then 4" . p? D s can cancel 
a propagator, whereas pf-2e) wn ^ &i ve r i se to an integral in 6 — 2e dimensions. It should 
be kept in mind that the specification given here is just a simple prescription relevant to 
practical calculations. The scheme is rigorously defined in Ref. [^7 . 

Awb,ioo P contains ultraviolet and infrared divergences. Ultraviolet divergences are re- 
moved after renormalization of the quark fields 

Ipbare = i^/Vren • (76) 

The renormalized amplitude is obtained as 

Awb,loop,ren = ( Z i/ ) -^Wb, loop, bare ■ (77) 



We have to renormalize the fields such that the residuum of the propagators is 1. For 
light quarks the appropriate renormalization constant is 1, due to a cancellation of UV- 
and IR- divergent parts. In more detail we have 



Z 4>,onshell,FD ~ 1 + 2 (AlT) 2 ^ F ^~^ UV + ^ IR ^ ' 

12 1 9 2 

Z iP,onshell,HV = 1 + o (A^\2 Cf (~^UV + W + &IR ~ Ur) , (75 



where A = 1/e — 7 + ln47r. Here we have indicated with a subscript UV or IR the origins 
of the divergent parts as well as the origin of additional finite terms which arise from a 
cancellation of 1/e-terms with terms of order e. For massive quarks we have 



Z i!,onshell,FD ~ 1 + o (A \2 ^ F [ ~3A UV - 5 + 3 In 



2(4tt) 2 1 V /i : 



2 



1 x t y^\^ i 



After renormalization of the quark field we have 

D Wb,loop,ren,FD ~ Wb,loop,bare,FD 2 I /I 2 / W/b ' born 

Bwb,loop,ren,FD = ^Wb,loop,bare,FD ■ (80) 
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The only divergences left in the renormalized amplitudes are of infrared origin. The infrared 
divergent parts are given by 



n Wb,loop 



m — sq7 i m 

iR \2~ ' 2' 2 '4 ~ JT 2 + 2 n 7i 



-2A ( -A + - 7 - - ln47r + - - In + - In — B WbMrn 



n Wb,loop 



IR = -2A^A + 7 -ln47r+^-ln-^^ W6)born . (81) 



We now turn our attention to Ward identities and finite renormalizations. In the 
't Hooft-Veltman scheme 75 anti-commutes with the first four Dirac matrices and commutes 
with the remaining ones. This treatment violates a Ward identity, which has to be restored 
by a finite renormalization. Explicitly, one splits the left-handed interaction into a vector 
(7 M ) and axial-vector \T ^ = 1/2(7^75 - 757^)] P art - 

1 1 1 

( 1_ 75) = ^ - 4 (7^75 - 757m) ■ (82) 

The Ward identity is restored by a finite renormalization of the axial-vector coupling: 

■pbare _ y -prenorm 
1 /i5 ~ axial, HVJin 1 ^5 

g 2 

Z axial, HV, fin = 1 + 4 7- ,~ Cp ■ (83) 

[4:71 ) 

There is some freedom in how we continue the left-handed coupling in the 't Hooft- 
Veltman scheme to D dimensions. For example, the expressions 

7 M ^( 1 -75), ^(1 + 75) 1^ ^(1 + 75)7m^( 1 -7s), ^7 M - ^ (7^75 - 757^) , (84) 

all agree in four dimensions, but differ in D dimensions. Of course this difference is 
compensated by the appropriate finite renormalization. Z axia ^ HV j in corresponds to the 
choice -1/4(7^75 - 7 5 7 M ). 

The four-dimensional scheme violates a Ward identity as well, which is restored by a 
finite renormalization of the left-handed coupling Ti e f t = 1/2^^(1 — 75), 

1 left - ^EW,FDfivX left 

g 2 

Z EW,FD,fin = 1 + x^ C f • (85) 

(47T y 



After the finite renormalization the renormalized amplitudes 

) 4 

t 4 



Awb,loop,ren,HV = \Z^,HVJ ^Wb,loop,HV 

^Wb,loop,ren,FD (4!fb) (ZEW,FD,fin) 2 ^-Wb,loop,FD , (86) 
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agree up to terms resulting from a different treatment of the infrared divergences. The 
axial-vector coupling in the 't Hooft-Veltman scheme has been renormalized according to 



-pbare 
1 /Lt5 



y -prenorm 
axial, HV, fin'- ^5 



In detail, we have for the finite parts of UV-origin in terms of g 2 / (47t) 2 CfB 



Born 



Bwb,loop,bare 

7-1 
coupling 

Zip 


FD 


-1 



HV vector 

-2 


+1 


H V axial 

+2 
-4 
+ 1 


Sum 


-1 


-1 


-1 



(88) 



The remaining differences are due to finite terms of infrared origin. Due to the universal 
structure of the infrared divergences we may relate the amplitudes calculated in the four- 
dimensional scheme to the ones in the 't Hooft-Veltman scheme 
are 

1 



8lB9fl- The relations 



B m 

D Wb,loop,HV 
D Wb,loop,HV 



B 



[1] 

Wb,loop,FD 



, D Wb,Born 



o[2] _ B [l) 

D Wb,loop,FD D Wb,Born ■ 



This completes our discussion on scheme-independence. In summary, we are able to 
perform the calculation in a four- dimensional scheme and to obtain from this result the 
amplitudes in the 't Hooft-Veltman scheme through simple and universal relations. 

We would like to comment on the original formulation of dimensional reduction. Di- 



mensional reduction |)0],|9T| differs from the four-dimensional scheme by use of the relation 



iW(4) 



P(D) ■ 1(4) , 



(90) 



which can be interpreted as "D is smaller than 4". The results in dimensional reduction 
can easily obtained from ours by dropping the terms Cq . The Cq 2e ''-terms yield finite 
terms related to UV-divergences. We note that the results in dimensional reduction cannot 
be related to ours, nor to the 't Hooft-Veltman scheme by a finite renormalization, since the 
term with Cq 2 ^ in B^ b loop is not proportional to Bwb,bom- Therefore the naive approach 
of dimensional reduction is not consistent. The situation may be cured at the expense 
of introducing additional scalar ghost particles f9"2|, however this spoils the calculational 
simplicity of the scheme. From a calculational point of view we prefer the scheme defined 
in Ref . |87| , since 6 — 2e-dimensional integrals are rather "inexpensive" to evaluate. 



4.3 Spin observables 

The helicity amplitudes contain the complete spin information for the single-top-quark 
processes. In the old-fashioned method, spin observables are calculated by inserting the 
spin projection operator 

1 

2 



u{p, s)u(p, s) 



(j) + m) - (1 + 75; 



(91) 
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into the matrix element squared. In Eq. ([91] ) s denotes a spin four- vector with s 2 = — 1, 
and p ■ s — 0. In the rest frame of the particle the spatial components of s point in the 
same direction as the spin of the particle |]9*3"|| . To make contact with this formulation, we 
first introduce the spin density matrix in the basis of Eq. (|53|) : 



P 



A(...,p+, ...)A%..,p+, ...) A(...,p+, ...)A*(...,p 
A(...,p~, ...)A%.,p+, ...) A(...,p~, ...)A*(...,p~ 



In addition we need the projection operator Eq. in the basis of Eq. 
obtained as 

-(q-\(jl> + m)(l-t)\q-) (q - \(J> + m)(l - 

-(<?+ |(/$ + m)(l-/%-> (g + + m)(l - 





(93) 



Spin observables are then calculated as 

Tr (P-p). (94) 

Note that the entries of the matrices are complex numbers, and that the spin vector s 
enters only through the matrix P. It is easily verified that this expression agrees with 
the one obtained from Eq. (|9l|). The spin summed result is recovered by replacing P in 
Eq. (|9~4"D with the unit matrix. 



5 Numerical results 

The inclusive NLO cross sections for s- and t-channel production of single-top-quarks were 



published in Refs. [[)6|] and |p8| , respectively. After we choose our numerical inputs, we 
compare to these older calculations, and update the inclusive cross sections using newer 
parton distribution functions (PDFs). 

In order to make a definite comparison to the older calculations, we reevaluate all 
results with the following parameters: For the mass of the W boson we use = 80.4 
GeV. For the top-quark mass we take m t = 175 GeV. In LO calculations we use CTEQ5L 



PDFs [§gj. In NLO calculations we use CTEQ5M1 PDFs with 2-loop running of a s . We 
define the electroweak coupling by g 2 = 8GpM^/ 'y/2, with Fermi coupling constant of 
Gp = 1.16639 x 10~ 5 GeV -2 . We consider pp collisions with center-of-momentum energy 
V~S = 1.8, 1.96, or 2.0 TeV (Tevatron), and pp collisions with a center-of-momentum energy 
^=14 TeV (LHC). 

Unlike the massive dipole formalism, the phase space slicing method of two cutoffs 
depends on explicit parameters 6 S and S c . In Figs. ^ and f| we show the s- and t-channel 
inclusive cross sections as a function of 5 S = 300 x 5 C . The logarithmic dependence cancels 
in the sum of the two- and three-body contributions and leaves terms proportional to 5 S 
and S c . By taking 5 S < a few xlO -3 , the cross sections converge to the updated analytic 
results. 
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Figure 3: The next-to- leading order single-top-quark total cross section for the s-channel at a 
y/S = 2 TeV proton-antiproton machine. The two- and three- body contributions, together with 
their sum, are shown as a function of the soft cutoff S s . The bottom enlargement shows the sum 
(open circles) relative to ±5% (dotted lines) of the analytic result (solid line). 
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Figure 4: The next-to-leading order single-top-quark total cross section for the t-channel at a 
y/S = 2 TeV proton-antiproton machine. The two- and three- body contributions, together with 
their sum, are shown as a function of the soft cutoff S s . The bottom enlargement shows the sum 
(open circles) relative to ±5% (dotted lines) of the analytic result (solid line). 
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Table 1: LO and NLO cross sections for single-top-quark production at the Tevatron and LHC 
for mt = 175 GeV. Cross sections are evaluated with CTEQ5L and CTEQ5M1 PDFs, and all 
scales set to mj. Errors include only Monte Carlo statistics. 



Process 






&NLO IPbJ 


s-channel 


l.o lev pp [t) 


U.ZOy 


n oon _i_n nno 
U.ooU ±U.UUz 




i.yo lev pp yij 








2 TeV pp (t) 


0.315 


0.463 ±0.002 




14 TeV pp (t) 


4.53 


6.55 ±0.03 




14 TeV pp (i) 


2.74 


4.07 ±0.02 


t-channel 


1.8 TeV pp (t) 


0.648 


0.702 ±0.003 




1.96 TeV pp (t) 


0.883 


0.959 ±0.002 




2 TeV pp (t) 


0.948 


1.029 ±0.004 




14 TeV pp (t) 


144.8 


152.6 ±0.6 




14 TeV pp (i) 


83.4 


90.0 ±0.5 



The LO and NLO cross sections for s-channel and t-channel single-top-quark production 
are listed in picobarns in Tables [l] and Since the Tevatron is a pp collider, the cross 
sections for top-quark (t) and antitop-quark (t) production are the same. The LHC is a 
pp collider, and hence the t and i cross sections are listed separately. Factorization (hf) 
and renormalization (/xr) scales are set equal. In Table [l] all scales are set to the top mass 
m t . In Table |] the s-channel cross sections are calculated using the invariant mass of the 
top-quark/antibottom-jet pair for the scale, /z 2 = M 4 | = (P t + Ph) 2 , where P t and are 
the four-momenta of the top quark and antibottom jet, respect ively.[] The t-channel cross 
section uses the double deep-inelastic-scattering (DDIS) scales, {if = Q 2 = ~(pb — Pt) 2 
for the light-quark line and [x\ = Q 2 + m 2 = —(p u — Pd) 2 + m 2 for the heavy-quark line, 
where Q 2 is the virtuality of the exchanged W boson (valid through NLO), and Pi are the 
four- momenta of the partons in Eq. (JI|). All LO cross sections are identical to those in 



Refs. [56] and [58] once updated to the CTEQ5 PDFs. 



We compare the NLO s-channel cross sections to a recoding of Refs. [9~5] and [Q. The 
results agree to within 1% for all scale choices. Of note is our use of M t i for the scale in 
Table [| rather than Q 2 , the virtuality of the W, which was used in Ref. |[56|| . These two 
scales are identical at LO and in initial-state corrections, but differ by the emitted gluon in 
final state corrections. Since we are interested in making cuts based on observables, such 
as M t i, we choose this as the scale. While the central value of any of the scale choices is 
very similar, the uncertainty is slightly larger at NLO using m t or M t j. In particular, if we 
vary the scales between M t i/2 and 2M t i (or m t /2 and 2m t ), we find the s-channel cross 
section varies by +7.8 - 6.9% at LO, and ±5.7 - 5.0% (±5.5 - 4.6%) at NLO. In contrast, 



using a scale of Q 2 , as in Ref. |56|, would predict a NLO scale uncertainty of ±4%. Given 



1 The choice of jet definition induces a cone-size dependence in the scale that always enters the cross 
section at one higher order in QCD. We have confirmed that this effect is numerically less than the overall 
scale uncertainty for any infrared-safe jet definition. 
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Table 2: LO and NLO cross sections for single-top-quark production at the Tevatron and LHC 
for nit = 175 GeV. Cross sections are evaluated with CTEQ5L and CTEQ5M1 PDFs, and M t i 
or the DDIS scales (m = Q 2 , \Xh = Q 2 + m 2 ), for s-channel or i-channel, respectively. Errors 
include only Monte Carlo statistics 



Process 


Vs 

V " 


oin (v>h) 

^ Li\J \ir KJ J 




s- channel 


1.8 TeV pp (t) 


0.244 


0.377 ±0.002 




1.96 TeV pp (t) 


0.287 


0.442 ±0.002 




2 TeV pp it) 


0.297 


0.459 ±0.002 




14 TeV pp (t) 


4.612 


6.56 ±0.03 




14 TeV pp (i) 


2.788 


4.09 ±0.02 


t-channel 


1.8 TeV pp (t) 


0.735 


0.725 ±0.003 




1.96 TeV pp (t) 


0.996 


0.990 ±0.002 




2 TeV pp (t) 


1.068 


1.062 ±0.004 




14 TeV pp (t) 


152.7 


155.9 ±0.6 




14 TeV pp (i) 


86.1 


90.7 ±0.5 



that we can probe a more restricted phase space with cuts, we take the conservative view 
that ±5.7% is an appropriate estimate of the NLO scale uncertainty at the Tevatron when 
looking at exclusive final states. At the LHC the scale uncertainty is less than ±2%. 

At LO and NLO, but not NNLO, color conservation forbids the exchange of a gluon 
between the light- and heavy-quark lines. Hence, the t-channel process may be factorized 
at NLO into two independent corrections that each resemble deep-inelastic-scattering. The 



NLO t-channel cross sections in Ref. were calculated using the scales suggested by this 
relation to double deep-inelastic-scattering. When using the DDIS scales, our results in 
Table match the updated Ref. [^HJ to better than 0.3% at the Tevatron, but are larger by 
2.9% at the LHC. The cross sections in Table [l] agree with the updated results of Ref. [58 



to within 1% when evaluated at fi = m t . In all cases, the dipole subtraction calculations 
and phase space slicing calculations agree within the statistical errors. 

A subtle issue arises in attempting to ascertain the effect of higher orders by varying 
the scale. In Ref. |58] only the scales in the vertex and PDFs of the heavy-quark line 
(b — > t) were varied because the corrections to the light-quark structure function are 
small.Q However, to the extent this process looks like double deep-inelastic-scattering, we 
expect a similarly small effect for the heavy-quark corrections as well. Indeed, the results 
for LO and NLO are nearly identical when using the DDIS scales, but differ significantly 
when using a fixed scale such as m t . In Table |3] we show the effects of varying the scales at 
the Tevatron (y/S = 2 TeV) together, and separately in the light- and heavy-quark lines. 
What we see is that varying the scales together at LO, particularly when using mt, severely 
underestimates the NLO correction. 



2 The NLO correction to the light-quark line effectively undoes the extraction of the NLO PDFs from 
the DIS data. 
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Table 3: Scale variation of the LO and NLO cross sections for i-channel single-top-quark pro- 
duction at the Tevatron (yS = 2 TeV). Variation in the light-quark and heavy-quark lines are 
listed as hi and Hhi respectively. Fixed (h = mt), and double deep- inelastic-scattering scales are 
shown separately. 







Hh & Hi 


Hh 


Hi 


LO t (mt) 
NLO t (mt) 


0.95 pb 
1.03 pb 


±1% 
±2.5% 


-7.5(57; 
+5.5/0 
-3.5q/ 
+4.0 /0 


+6.7q/ 
-5.8/0 

±1% 


LO t (DDIS) 
NLO t (DDIS) 


1.07 pb 
1.06 pb 


t° 2 A % 

±3.5% 


-7.2q/ 
+5.2/0 

+1% 


+8 0/ 
-6.8 /0 

±0.6% 



The reason for the underestimate in varying the scales together is a series of accidental 
cancellations that are driven by the range of proton momentum fraction probed at the 
Tevatron. For a top mass of 175 GeV and a machine of around 2 TeV, the typical x ~ 0.1. 
For scales also around 100-200 GeV, the h PDF happens to increase with increasing scale, 
whereas the valence quarks decrease. This may be seen in the opposite signs in the last 
two columns of Table |^. The net effect is that to estimate the uncertainty, we must 
vary the scales independently, and then add them in quadrature.^ This leads to a LO 
uncertainty of ~ ±10%, and a NLO uncertainty for the Tevatron of ±4%. At the LHC, 
the NLO uncertainty is ±3%. These uncertainties are consistent with expected higher-order 
corrections, and with both fixed and DDIS scale choices. 

While the s-channel cross section has only changed by a couple of percent from p6| , 
the t-channel cross section is 13% smaller at the Tevatron than appears in []58[| . The shift 
in the NLO i-channel cross section is due to the correction of bugs that appeared in the 
evolution of the gluon PDF in all of the older MRS and CTEQ NLO PDFs |9|. Because 



the b PDF is constructed almost entirely out of gluon evolution [^,|97|], the effect of the 
bug is greatly enhanced in all processes where there is a c or b quark in the initial state. 
The programming bug is corrected in MRS99 (updated) |43| and CTEQ5M1 p2| . 

The central goal of our calculations is not to recalculate inclusive cross sections, but 
to provide full momentum and spin dependent distributions with the option for arbitrary 
cuts. Detailed analyses of the phenomenological issues concerning these distributions will 
be presented elsewhere |98| . Here we restrict ourselves to a comparison of the results using 



the phase space slicing method and massive dipole formalism. 

For our comparisons we reconstruct jets using a kx cluster algorithm [99] with AR = 1. 
We define a simple detector acceptance by assuming that only jets with p T > 20 GeV and 
1 77 1 < 2 are observed. We calculate all cross sections at the scale h = m t = 175 GeV, and 
at a 2 TeV pp collider. 

In Figs. [5] and |] we present the NLO transverse momentum pxb and pseudorapidity 
r)b distributions of the 6-jet in s-channel production of a top quark, with cuts based on 

3 The light- and heavy-quark corrections are factorizable through NLO, and hence are only weakly 
correlated through the evolution of the PDFs. 
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the "jet veto" search strategy in Ref. ||. We accept events where only the top-quark and 
6-jet pass the cuts above, and any additional jets are either soft (pTj < 20 GeV) or are 
outside the simple detector (\r]j\ > 2). The phase space slicing method and massive dipole 
formalism give identical distributions, even in this region of restricted phase space. 




50 100 150 200 

PTb (GeV) 



Figure 5: Transverse momentum distribution of the 6-jet at NLO in s-channel production of a 
top quark at the Tevatron (\/ r S = 2 TeV) after cuts. Phase space slicing (PSS) results (solid) 
and massive dipole formalism (MDF) results (dashed) are both shown. 

In Figs. [7] and |8| we present the NLO transverse momentum pxj 1 and pseudorapidity 
rjj-L distributions of the highest-pr jet in t-channel production of a top quark, using the 
same "jet veto" search strategy as above. We accept events where only the top-quark and 
highest-pr jet pass the cuts above, and any additional jets are either soft ij>Tj 2 < 20 GeV) 
or are outside the simple detector (|r/ j2 | > 2). The phase space slicing method and massive 
dipole formalism predict identical distributions at NLO. 

6 Conclusions 

We present three independent calculations of the fully differential production of a single 
top quark plus one jet at next-to-leading order in hadronic collisions. At this order in 
QCD the cross sections factorize into two non-interfering production modes that may be 
identified by the s-channel or t-channel exchange of a IF boson. The s-channel cross 
section is characterized by having a typical final state of a top quark and a bottom-quark 
jet. The jet in the t-channel cross section tends to be somewhat more forward, and rarely 
contains a bottom quark. New physics scenarios tend to effect these production modes 
differently, and hence an accurate measurement and theoretical calculation of both modes 
is desirable. Further, the CKM matrix element \Vtb\ 2 is directly proportional to the cross 
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Figure 6: Pseudorapidity distribution of the 6-jet at NLO in s-channel production of a top quark 
at the Tevatron (y/S = 2 TeV) after cuts. Phase space slicing (PSS) results (solid) and massive 
dipole formalism (MDF) results (dashed) are both shown. 



section, and the ability to extract V t b is ultimately limited by our ability to predict the 
measured exclusive cross sections with experimental cuts. 

The total cross sections for s- and t-channel production are updated with CTEQ5 
parton distribution functions for runs I and II of the Tevatron, and for the LHC The 
t-channel cross section is smaller than previously published by 13%. This decrease is due 
entirely to the correction of a bug in the parton distribution functions. If the Tevatron 
continues to run at 1.96 TeV instead of 2 TeV, the s(t)-channel production modes will 
produce 4% (7%) fewer events than previously expected. We discuss some subtleties in 
estimating the effects of higher-order corrections on these cross sections. 

We show that the phase space slicing methods of one and two cutoffs and the massive 
dipole formalism produce the same jet distributions and cross sections. The dipole cal- 
culation retains the full spin dependence of the external particles, and so may be used to 
predict spin-dependent correlations of the partons. We discuss elsewhere |)8| the effects of 
top-quark mass, scales, jet definitions, and parton distribution functions on the shapes of 
the NLO distributions. 
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Figure 7: Transverse momentum distribution of the highest-py jet at NLO in t-channel 
production of a top quark at the Tevatron (y/S = 2 TeV) after cuts. Phase space slicing (PSS) 
results (solid) and massive dipole formalism (MDF) results (dashed) are both shown. 
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A Scalar integrals 

We calculate the integrals in D = 4 — 2e dimensions and use the notation 

A = i-7 + ln4vr. (95) 

£ 

All integrals are calculated in the Euclidean region (invariants p 2 < and masses m 2 > 0). 
These expressions can then be continued analytically to the regions of interest, using the 
substitution — s — > — s — ie (where e denotes a small parameter, not to be confused with 
the one appearing in dimensional regularization) and the formula 

\n(-s-ie) = \n(\s\) - me(s) . (96) 

Massive integrals: The tadpole is given for m 2 > by 

Mm >) = ( 4,)V- D /(S ? ra--'" 2r (- 1+£ )( 4 ^(S" 1 

= m 2 U + l-ln 7 ^) +0(e). (97) 
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Figure 8: Pseudorapidity distribution of the highest-pr jet (ji) at NLO in i-channel production 
of a top quark at the Tevatron (y/S = 2 TeV) after cuts. Phase space slicing (PSS) results (solid) 
and massive dipole formalism (MDF) results (dashed) are both shown. 



The bubble with one internal mass is given for — p 2 > 0, m 2 > by 

d D k 1 



4V,m 2 ) = (4tt)V- D / 



(27r) D i k 2 ((k - p) 2 - m 2 ) 



= A + 2 -In In— + 0(e) . (98) 

fi z —p A m z — p z 

On-shell integrals: The bubble with one internal mass is given for p 2 = m 2 > by 

B^im 2 ) = (4ir) 2 ^- D [ — L 

0V ; V ;P J (2vr) D i k 2 ((k-p) 2 - m 2 ) 

9 

r ji * 

= A + 2-ln— +0(e). (99) 
A* 

The triangle with two external masses and one internal mass, between the two external 
massive lines is given for p\ = 0, (pi + P2) 2 = m 2 > by 

Ci" } (A, m 2 ) = (to) / {2 ^ Di kl(k _ pi y ({k _ pi _ p2 y_ m , ) 

- M^{-i4('"H-rf)4' n «™ 2 ») 



+1 W(^) - i ( ro » - ^) + Li 2 (-j^) } + OW 
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1 a 2 a ( 1 1 1 . 1 . m 2 . m 2 -pV\ 1 , , , x2 
A + A --7 + -ln4vr + - In — + In - - (7 - ln47r) 



m 2 - p\ \ 2 V 2 ' 2 2 /z 2 m 2 / 4 

^(5)-H^)M^)-s} +ow - (ioo) 

Integrals with no internal masses: The bubble with — p 2 > is given by 

2 

= A + 2-ln^- + 0(e). (101) 

The triangle with one external mass p\ = p\ = 0, — s = —2pip 2 > is given by 

d D k 1 



^ 6) ( S ) = (4 7 r)V- D / 



(27r) D i k 2 (k-p 1 ) 2 (k-p 1 -p 2 ) 



2 



n rfl + «)!*(!-«) 
s£ 2 r(i-2s) 1 ' l^ 2 / 1 

i (a' + a ( 7 - 1„ 4. - m =£) + 1 (7 - m4,)> + 1 fe) - ^ ] + 0(e) 



(102) 



Six-dimensional integrals: The triangle with an additional power of k 2 _ 2e in the numer- 
ator is equivalent to an integral in D = 6 — 2e dimensions and yields for the massless and 
massive case 



Ct 2£) = (4vt)V- D / d ° k kl 



(~2e) 



2^ 



(2n) D i k 2 (k — Pi) 2 ((k — p\— P2) 2 — m 

= -\ + 0{e). (103) 
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